Shape Coexistence in 186 Pb: Beyond-mean-field description 
by configuration mixing of symmetry restored wave functions 
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We study shape coexistence in 186 Pb using configuration mixing of angular-momentum and 
particle-number projected self-consistent mean-field states. The same Skyrme interaction SLy6 
is used everywhere in connection with a density- dependent, zero-range pairing force. The model 
predicts coexisting spherical, prolate and oblate + states at low energy. 
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The behavior of shell effects away from the valley of 
stability is a topic of very active investigation, both the- 
oretically and experimentally For light nuclei, the spher- 
ical N — 20 and N = 28 shells disappear in neutron-rich 
isotopes, leading to strongly deformed ground states and 
large B(E2) transition probabilities between the first 2 + 
state and the ground state In contrast, the magic 
proton number Z = 82 is particularly strong and its in- 
fluence persists even in very neutron-deficient nuclei. The 
ground state of Pb isotopes is known to be spherical down 
to 182 Pb @. The weakening of the magicity of the Z = 82 
shell manifests itself through the appearance of low-lying 
+ states 3J. At least one low-lying, excited + level 
has been observed in all even-even Pb isotopes between 
A = 182 and 194 at excitation energies below 1 MeV, the 
most extreme cases being 188 Pb and 186 Pb 0,0 with two 
excited + states below 700 keV. 

Two different kinds of models have been invoked to 
explain the coexistence of several + states at low en- 
ergy . In a shell model picture , the first excited + 
level observed from 202 Pb down to 186 Pb is interpreted as 
a two-quasiparticle proton configuration (7r/i9/2) 2 , while 
the second one in 188 Pb and 186 Pb as well as the first 
+ state in 184 Pb are understood as a four-quasiparticle 
configuration (7r/i9/2) 4 . In this picture, neutrons and 
protons outside the inert core interact through pairing 
and quadrupole interactions to generate deformed struc- 
tures. Such a model requires a drastic truncation of the 
configuration space. Up to now, it has only been applied 
in rather schematic and qualitative ways. 

In mean-field models, the + states observed at low 
energies are associated with coexisting energy minima 
which appear for different values of the axial quadrupole 
moment jjj . The ground state corresponds to the spheri- 
cal minimum and the excited + level to a deformed state 
with an oblate (in the heaviest Pb isotopes) or a prolate 
(in 184 Pb up to 188 Pb) shape. 

However, shape coexistence in the neutron-deficient 
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Pb region cannot be described on the level of mean-field 
models in a fully satisfactory way. The minima obtained 
as a function of the quadrupole moment are rather shal- 
low and dynamical effects such as quadrupole vibrations 
mayaffect the very existence of these minima. Tajima et 
al. Q and, more recently, Chasman et al. have studied 
the quadrupole dynamics of Pb isotopes by performing 
a configuration mixing of mean-field states with different 
axial quadrupole moments. Their results support the 
interpretation of the excited + states as deformed min- 
ima. The lowest excited levels obtained in the configura- 
tion mixing calculation have average deformations close 
to that of the mean-field minima. However, the calcu- 
lated excitation energies overestimate the experimental 
values. Diabatic effects have been studied by Tajima 
et al. who have included, for each axial quadrupole mo- 
ment, the lowest Hartree-Fock+BCS (HFBCS) configu- 
ration and the two-quasiparticle deformed proton config- 
urations (7r/i9/2) 2 . Tajima et al. have shown that these 
configurations do not influence the configuration-mixing 
results significantly, and that they can be neglected. 

The experimental data on neutron deficient Pb iso- 
topes are not limited to a few + states. Rotational bands 
have also been observed whose properties have served to 
interpret the excited + state as associated with oblate 
and prolate deformations. Transition probabilities be- 
tween the levels are also known in some cases. It seems, 
therefore, highly desirable to apply the configuration- 
mixing method that we have recently developed jlCj to 
Pb isotopes. This method treats simultaneously the most 
important symmetry restorations and the mixing with 
respect to a collective variable. Here, we present an ap- 
plication to 186 Pb. This isotope has the unique property 
of having + levels as its lowest three states, with the 
excitation energy of the second and third + also being 
the lowest among the known Pb isotopes While the 
ground state is assumed to be spherical, the + states 
observed at 532 keV and 650 keV are interpreted as cor- 
responding to oblate and prolate configurations. 

The "projected" configuration mixing of mean-field 
wave functions performed here has several goals. The 
particle-number projection removes unwanted contribu- 
tions coming from states with different particle numbers, 
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which are an artifact of the BCS approach. The angular 
momentum projection separates the contribution from 
different angular momenta to the mean-field states and 
generates wave functions in the laboratory frame with 
good angular momentum. Finally, the variational con- 
figuration mixing with respect to a collective coordinate, 
the axial quadrupole moment in this work, removes the 
contributions to the ground state coming from collective 
vibrations, and simultaneously provides the excitation 
spectrum corresponding to this mode. 

The starting point of our method is a set of indepen- 
dent HFBCS wave functions \q) generated by mean-field 
calculations with a constraint on a collective coordinate 
q. Such mean-field states break several symmetries of the 
exact many-body states. Wave functions with good an- 
gular momentum and particle numbers are obtained by 
the restoration of rotational and particle- number symme- 
try on \q): 

l JM «) = JfT,9K P MKPzPN\q), (1) 

where N is a normalization factor. P'l IK , Pn, Pz are pro- 
jectors onto the angular momentum J with projection M 
along the laboratory z-axis, neutron number TV and pro- 
ton number Z, respectively. We impose axial symmetry 
and time reversal invariance and, therefore, K can only 
be and we shall omit the coefficient g J K = Sko- This 
prescription excludes the description of 7 bands where 
K = 2. 

A variational configuration mixing on the collective 
variable q is then performed for each angular momentum 

\JMk)=J2fk M (<l)\JMq}. (2) 

Q 

The weight functions // (<?) are determined by requiring 
that the expectation value of the energy 

JM = (JMk\H\JMk) 

k (JMk\JMk) [ ' 

is stationary with respect to an arbitrary variation 
8f£ M (q). This prescription leads to the discretized Hill- 
Wheeler equation . Such a secular problem amounts 
to a restricted variation after projection in the set of 
states obtained for different values of the collective vari- 
able q. Collective wave functions in the basis of the in- 
trinsic states are then obtained from the set of weight 
functions f(! M (q) by a basis transformation In 
\JMk), the weight of each mean-field state \q) is given 
by: 

g( M (q) = (JMk\q), (4) 

Since the collective states | JMk) have good angular mo- 
mentum, quadrupole moments and transition probabili- 
ties can be determined directly in the laboratory frame 
of reference without further approximations. 
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FIG. 1: Particle-number projected ("mean field") and 
particle-number and angular-momentum projected potential 
energy curves up to J = 10 for 186 Pb as a function of the mass 
quadrupole moment in barn (upper axis) or, equivalently, in 
terms of (h (lower axis). The energy reference is that of the 
projected spherical mean-field state. 



The same effective interaction is used to generate the 
mean-field wave functions and to perform the configu- 
ration mixing calculation. We have chosen the Skyrme 
interaction SLy6 in the mean-field channel |l2^ and a 
density-dependent, zero-range force as defined in |l3j |. in 
the pairing channel. The pairing equations are solved us- 
ing the Lipkin-Nogami prescription, as done in |Toj . The 
two-body center-of-mass correction is self-consistently 
included in the interaction SLy6 |12| . However, in 
the present calculations, it is included a posteriori at 
the mean-field level as well as in the projection and 
configuration- mixing calculations. 

In figure H the deformation energy of 186 Pb is plotted 
before and after projection on angular momentum. All 
curves are drawn versus the intrinsic axial quadrupole 
moment of the unprojected mean-field states. As pro- 
jected J = states are spherical, this "quadrupole mo- 
ment" is only a convenient way to label the projected 
states. The curve labeled "mean-field" plots the defor- 
mation energy after particle-number projection only. It 
exhibits a spherical global minimum as well as local min- 
ima at prolate and oblate deformations. While the de- 
formation energy of the prolate minimum fortuitously re- 
produces the experimental value of 0.650 MeV for the 
prolate + state, the 1.1 MeV deformation energy of the 
oblate minimum overestimates the experimental value of 
0.532 MeV for the oblate 0+ level. A fourth, very shal- 
low, minimum can be seen at a deformation fa ~ 0.5; 
it is too shallow to be safely associated with a physical 
state. 

The energy curves obtained after angular momentum 
projection are also shown in figure ^ At moderate de- 
formations, around the prolate and oblate minima, the 
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FIG. 2: Spectrum of the lowest positive parity bands with 
even angular momentum and K — 0, as a function of the de- 
formation (see text). The angular momentum projected en- 
ergy curves are shown for comparison. The energy reference 
is that of the calculated 0^ ground state. 



mean-field states are dominated by angular momentum 
components with J > 8. This is reflected in the fact that 
all projected energy curves are far below the mean-field 
one. The spherical mean-field state is rotationally in- 
variant and, therefore, contributes to J = only. Two 
minima appear at small deformations, around 02 = ±0.1. 
They do not correspond to two different states, but to the 
correlated spherical state (see below). For larger prolate 
and oblate deformations, the energy difference between 
the mean-field and J = curves stays nearly constant. 
The prolate and oblate mean-field minima are present in 
all the projected energy curves. Angular momentum pro- 
jection reduces the energy difference between the spher- 
ical (|/32 1 ~ 0.1) and deformed minima to 0.2 MeV for 
the prolate and 0.68 MeV for the oblate well. While 
the prolate potential well is pronounced for all angular 
momenta, the oblate one now becomes very shallow for 
J = 0. 

The excitation energies E% M of the collective states 
\JMk) obtained from the configuration mixing calcula- 
tion are shown in figure [3 Each of these states is repre- 
sented by a horizontal bar drawn at the average intrinsic 
deformation J^g/M?) lfffe M ('?)| 2 ; where ,$2(9) is the de- 
formation of the mean-field state. The excitation spec- 
trum is divided into bands associated with different de- 
formations. Configuration mixing lowers the energy of 
the lowest collective states with respect to the projected 
energy curves. The energy gain is the largest for the 
ground state, hence increasing the excitation energies of 
the prolate and oblate + levels. 

The corresponding collective wave functions g£ M (q) 
are presented in figure Their square gives the weight 
of each mean-field state \q) in the collective state | JMk). 

The ground state wave function is spread in a similar 
way on both oblate and prolate sides with a zero value for 



the average P2 deformation. The wave functions of the 
first two excited + states are strongly peaked at either 
prolate (0^~) or oblate deformation (O3 ), with their tails 
extending into the spherical well. For higher J values, 
the shape of the wave functions confirms their assign- 
ment to oblate and prolate bands, as already hinted so 
in figure|21 As there is no spherical well for J > states, 
their wave functions mix only prolate and oblate configu- 
rations. Starting with J = 4, all levels are strongly local- 
ized and are predominantly either prolate or oblate. The 
shapes of the O4", 2^", and 4^ wave functions suggest their 
interpretation as a rotational band built on a [3 vibration 
within the prolate well, while the wave function of the l\ 
state indicates that it corresponds to a vibrational state, 
which is spread over the entire potential well. 

The calculated spectrum is compared with the exper- 
imental data in figure The excitation energy of the 
prolate + state at 0.55 MeV is very close to the ex- 
perimental value. In contrast, the excitation energy of 
the oblate + level is largely overestimated. The experi- 
mental data even suggest that the oblate state is slightly 
below the prolate one. A nice result from the calculations 
is that the structure of the first three + levels is dom- 
inated by spherical, prolate and oblate configurations, 
respectively, and this supports the interpretation of the 
experimental data in terms of shape coexistence. This 
feature could not have been guessed from the deforma- 
tion energy curves (see figure ^) , where the oblate well 
has a depth of only 500 keV in contrast with the prolate 
one of 1 MeV. The excitation energies of the first two 
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FIG. 3: GCM wave functions of the lowest | Jk) states. Solid 
lines denote spherical, long-dashed lines oblate, dashed lines 
prolate, and dotted lines the (3 band in the prolate well. 
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FIG. 4: Comparison between the calculated excitation ener- 
gies and the available experimental data for low-lying states 
in 186 Pb. From the left to the right the spectra show oblate, 
spherical and prolate bands. 



excited + states are even quite close to the energy dif- 
ferences between the deformed minima and the spherical 
minimum of the mean-field deformation energy curve. 

Both experimentally and theoretically, all bands ex- 
hibit a rotational behavior, with the exception of the 
E a + — E 2 + energy difference which is too small. This 
can be understood from the stronger state mixing for the 
J = than for higher J values which is observed in the 
calculations. For the prolate band, however, the displace- 
ment from a rotational behavior remains too small. This 
is probably a consequence of the overestimated energy of 
the oblate band head which reduces the mixing between 
the deformed configurations. 

Calculated transition probabilities for J > 2 states 
confirm the separation of the excited states into rota- 
tional bands with very small B(E2) transitions between 
them. While the transition quadrupole moments, Qo, of 
the oblate (Qo ~ —600 e fm 2 or /?2 ~ —0.2) and prolate 
(Qo ~ 1000 e fm 2 or 02 « 0.34) bands slowly grow with 
angular momentum, the deformation of the third rota- 
tional band stays nearly constant at about Qo ~ 1400 e 
fm 2 (02 ~ 0.49), in agreement with the systematics of the 
minima in the projected energy curves of figure ^ The 
B(E2) values for the in and out of band 2 + — > + tran- 
sitions confirm that the low-lying + states are indeed 
mixed. 

Our results strongly support the interpretation of the 
Pb isotopes spectra as evidence for shape coexistence. 
There remains, however, a significant over-estimation of 



one of the band's excitation energy. This could be due 
to several ingredients of the model: 

• the effective mean-field interaction; small differences 
between interactions (surface tension, spin-orbit strength 
. . . ) shift the relative energies of the various coexisting 
minima at the mean-field level 0, 0] . In a calculation 
with the Skyrmc SLy4 interaction, the prolate and oblate 
0+ states are pushed up to 1.05 MeV and 1.39 MeV, 
respectively, as can be expected from the overall stiffer 
energy surface of this interaction |T(| . 

• the strength and the form factor of the pairing in- 
teraction; a test with a reduced pairing strength (—1100 
MeV fm 3 ) shows that the energies of the prolate and 
oblate minima of the deformation energy curves are re- 
duced to 0.2 MeV and 0.65 MeV, respectively. 

• the configuration space used in the configuration 
mixing; to test this possible source of error, we have 
enlarged the space by including the oblate (7r/i9/2) 2 
two-quasiparticle proton configurations, as was done by 
Tajima et al. [8j. As in this work, the results are changed 
by at most 100 keV. 

• the inclusion of triaxial quadrupole configurations; 
projection on J becomes much heavier numerically, and 
this is still beyond present numerical possibilities. 

• generalized interaction for calculations beyond mean 
field; most mean-field interactions depend on the one- 
body density. It is known since the 70s that this den- 
sity dependence can have two different origins: either a 
three-body force or a resummation of (short-range) cor- 
relations. To generalize this dependence for the non- 
diagonal matrix elements appearing beyond the mean- 
field approximation, we have chosen the generalisation 
stemming from a three-body interaction. Resummation 
of correlations beyond mean-field gives rise to another 
generalisation of the Skyrme force [l7|. The study of 
shape coexistence in nuclei like the Pb isotopes could be 
a good place to determine the merits of both generalisa- 
tions of the Skyrme force. 
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